o 
o 
o 

+-> 

q ; On the motion of rotating bodies 

in field gravity theory and general relativity 

7— I . 
> 

in , 

q , Yurij Baryshev 

O ' 

Institute of Astronomy, St. Petersburg State University, 

o : 

O ■ Staryj Peterhoff, 198504 St. Petersburg, Russia, 

qh e-mail: yuba@YB3255.spb.edu 



X 



Abstract 

On the basis of Lagrangian formalism of relativistic field theory post-Newtonian 
equations of motion for a rotating body are derived in the frame of Fcynman's 
quantum field gravity theory (FGT) and compared with corresponding geodesic 
equations in general relativity (GR). It is shown that in FGT the trajectory of 
a rotating test body does not depend on a choice of a coordinate system. The 
equation of translational motion of a gyroscope is applied to description of laboratory 
experiments with free falling rotating bodies and rotating bodies on a balance scale. 
Post-Newtonian relativistic effect of periodical modulation of the orbital motion of a 
rotating body is discussed for the case of planets of the solar system and for binary 
pulsars PSR B1913+16 and PSR B1259-63. In the case of binary pulsars with known 
spin orientations this effect gives a possibility to measure radiuses of neutron stars. 
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1 Introduction 



The discovery and continuous observations of binary pulsars has opened a new possibil- 
ities for testing relativistic gravity theories |ffl]. Timing observations of a pulsar in orbit 
with a stellar-mass companion require practical application of different relativistic grav- 
ity effects such as advance of periastron, gravitational redshift, time dilation, bending 
of light, Lense-Thirring precession, geodetic precession, and also gravitational radiation 
(see e.g. § - §). 

All this effects are usually calculated in the frame of general relativity (GR) which 
is the geometrical approach to gravity. Observations of binary pulsars verified GR with 
high accuracy, however it is still important to apply other gravity theories for binary 
pulsars observations |5|, ||. 

Here it is considered the field gravity theory which elements have been formulated 



in 60-th by Feynman |T(|, Thirring |Hj |, Kalman [12] and others (for recent discussion 
see Straumann jnj, Baryshev |H|] ). In Lectures on Gravitation |l(]] Feynman gave 
a quantum relativistic field description of gravity similar to other fundamental forces 
(strong, weak, electromagnetic). Feynman's field gravity theory (FGT) is based on 
the quantum field theory and presents an alternative non- geometrical understanding the 
physics of gravitational interaction. He emphasized that though the theory of gravitation 
"has both a field and a geometrical interpretation" , "the geometrical interpretation is not 
really necessary or essential to physics" (p0|], p. 113). According to Feynman, the gravity 
force between two masses is caused by the exchange of gravitons which are mediators of 
the gravitational interaction and actually represent the quanta of the relativistic tensor 
field i)) lk in Minkowski space rf k . 

The fact that the field equations in FGT and GR is the same in linear approximation 
and that in FGT there is an iteration procedure which saves gauge-invariance and leads to 
the same non- linear field equations, is usually interpreted as complete identity of GR and 
FGT. However, after the iteration is fulfilled and geometrical interpretation is accepted, 
the new puzzling conceptual problems appear, among them the spacetime horizons and 
the break down of the energy-momentum conservation in geometrical approach (see 
review by Straumann |13|). Moreover the iteration procedure in the frame of FGT is 
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not uniquely denned operation because it requires to fix an expression for the energy 
momentum tensor of the gravitational field which is not uniquely defined in Lagrangian 
formalism and needs additional physical assumptions such as symmetry, positiveness of 
the field energy, traceless energy-momentum tensor for massless fields (see Landau & 
Lifshitz ||, Bogolubov & Shirkov @, Baryshev @). 



All classical PN relativistic gravity effects are the same in FGT and GR (see e.g. 
[11, 14, Eo|). However here it is shown that there is a post-Newtonian effect which 



demonstrate a difference between FGT and GR and which may be tested by timing ob- 
servations of binary pulsars and by future laboratory experiments with rotating masses, 
e.g. measuring a balance between differently rotating bodies. This post-Newtonian 
effect leads to periodic modulation of the orbital motion of a pulsar, caused by the 
velocity dependence of gravity force acting on a rotating body. In GR this effect is a 
coordinate-dependent concept and hence could not be uniquely defined, while in FGT 
there is definite coordinate-independent value for this relativistic gravity phenomenon. 



The effect was firstly discussed in [17] as a prediction of the field gravity theory for 
the Stanford gyroscope experiment [18]. Here this effect is discussed in connection with 
binary pulsars timing observations and formulation of the Galileo-like experiments with 
free falling rotating bodies or measuring a balance between rotating bodies. 

In sec. 2 post-Newtonian equations of motion of relativistic test particles in relativis- 
tic tensor gravity field are discussed and compared with geodesic equations of general 
relativity. In sec. 3 equations of motion for a gyroscope orbiting around a central mass 
are derived and applied to the cases of laboratory experiments, planets of the solar sys- 
tem and binary pulsars as tests of the velocity dependence of the gravity force acting on 
rotating bodies. Conclusions are presented in sec. 4. 

2 Post-Newtonian equations of motion for test particles in 
field gravity theory and general relativity 



3 



2.1 Kalman's equation of motion in FGT 

The basic principles of FGT are the same as for other relativistic quantum fields. These 
include the Minkowski space, the fundamental role of the inertial frames, the concepts 
of force and energy-momentum tensor, and also quantum uncertainty principle and the 
many-path approach. 

Following (nj let us consider the motion of relativistic test particle with rest-mass 
mo, 4- velocity u l , 3- velocity v in the gravitational field described by the symmetric tensor 
potential tp on flat Minkowski space-time where the Cartesian coordinates always exist 
and the metric tensor rf k = diag(l, —1, —1, —1) (here we accept notations of the text- 
book pp. 

To derive the equation of motion in FGT one may use the stationary action principle 
in the form of the sum of the free particles and the interaction parts: 

8S = J (A {p) + A (mt) )dn) = (1) 

where d£l is the element of 4- volume and the variation of the action is made with respect 
to the particle trajectories 5x l for fixed gravitational potential ijj lk . 
The free particle Lagrangian is 

A w = -rfofcTg (2) 

and the interaction Lagrangian is 

A(int) = --^TpikTfy (3) 

where the energy-momentum tensor (EMT) of the point particle is 

2$ = m c 2 5(r- ? P ){1 - ^} V Vu fe (4) 

Inserting equations (||),(||) into equation(|l]) and taking into account that ds 2 = dxidx 1 
we get 

J (rriQc5{\J dxidx 1 ) + ~^ ik ~j d =j f ) = ( 5 ) 
Opening the variation and integrating by parts, and taking into account that the varia- 
tion is made for the fixed values of the integration limits, we find 

/ m cdu i 5dx i + ^^-d(u k ^i k )Sx i - —d(ipi k u l u k u i )5x i - —vPS^u-dx 1 = (6) 
J c c c 



Consider also that 

du 1 

du 1 = ——ds; dx l = u l ds; 5ip k i = "4>kil& x ■> 
ds 

d(u l u h Uiifjik) = u l UiU k 7pi ktn dx n + u l u k ipi k dui + 2tpi k u l l Uidu k 

finally we get the following equation of motion for test particles in the field gravity theory 
[19]: 

4^ = -m Bl lU k u l (7) 
where p k = rriQCu k is the 4-momentum of the test particle, = and 

1 2 2 

A\ = (1 - ^ nl u n u l )rf k - -^i) kn u n u l + -^f k (8) 

2 11 

B kl = ^k,l ~ ^kl ~ ^kl,nU n U l (9) 

The rest mass of the test particle can be canceled, hence in the field gravity theory 
f^in = mg = mo without initial equivalence postulate as was emphasized by Thirring 
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The equation (0) is identical to the equation of motion derived by Kalman |l2| in 
another way, by considering the relativistic Lagrange function L defined as S = J L^- 
and relativistic Euler equation: 

d . dL u , dL dL . . 

Ts^ u - L)ui + M=-W (10) 

Inserting in equation ( JToj ) the following expression of the relativistic Lagrange function 

L = -m c -m ipik-: — -j- (11) 
ds ds 



one gets Kalman's equations of motion [12], which may be transformed into equation 
©■ 

2.2 PN equations of motion in FGT 

In post-Newtonian approximation we keep terms an order of v 2 /c 2 ~ ^pn/c 2 <C 1 in 
equation (0). For PN accuracy we need calculations of the ijj 00 component with the 
same order, while other components of the tensor gravitational potential tp can be 
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calculated in Newtonian approximation. Under these assumptions from equation (||) for 
(i = a) we get the expression for the PN 3-dimensional gravity force (which we shall call 
the Poincare gravity force remembering his pioneer work in 1905 on relativistic gravity 
force in flat space-time): 



" / ' 3r 4> -ft, w , v -± . - 

Fpoincare = 37 = -m {(l + --=f + 3^)V0 - 3-(- • V0)| 

at 2 c z c z c c 



-™ {3-^- - 2-^| + 2(- x roi*)} (12) 
c c at cm c 

where ^ = V> 00 , * = ip 0a = "V'Oa- 

An important advance of FGT is that in Minkowski space according to Noether's 
theorem there is well-defined concept of the energy-momentum tensor and its conser- 
vation in the form which is usual for relativistic quantum fields. Hence the energy of 
gravity field is localized positive physical quantity and is given by the 00-component of 
the energy- momentum tensor (and not pseudo-tensor as in GR). For post-Newtonian 
applications it is sufficient to use the weak field approximation for the calculation of 
the gravity field EMT, which gives in the case of a static spherically symmetric massive 
body the following value for the 00-component (]ll|, |20tl): 

^=^=dgcw 3 (13) 

where ipjy is Newtonian potential. The energy of the gravity field gives the following PN 



correction to the 00-component of tensor gravitational potential ([11, p0| ) : 

<p = i; ™ = VN + l?k (14) 



Hence in the frame of FGT the energy of the gravity field is observable quantity which 
may be measured by observations of test particles motion. Take into account the ex- 
pression ( |l4]) for the 00-component of the gravitational potential, we get corresponding 
PN 3-acceleration of the test particle : 



^ = -< 1 + ^ +4 ^ V ^ + 4 c ( c- V 



+ '-cHt- 2 78i +2 ^^ (15) 
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^From the (i = 0) component of equation (^) it follows the expression for the work 
of the Poincare force: 

dEk -j? - in m**^ d# n^ 

~^=v- F Pancare = -m v ■ {(1 - + 3- 2 W - 3-— + 2—} (16) 

An important particular case is the static spherically symmetric weak gravitational 
field for which ^ = 0, d(j)/dt = 0, ijj lk = diag((j),ipN ,tfN jVat) hence PN 3-acceleration 
will have the simple form: 



{~tt)fgt = - 1 + -J + 4^)V(p N + 4- - ■ V^tv 17 



^From equation of motion ( [i7[ ) it is clear that the acceleration of a test particle de- 
pends on the value and the direction of its velocity and this is coordinate independent 
relativistic gravity effect. 

For circular motion v _L V</3jv hence PN 3-acceleration is 



(fw = -(l + J + 4^)V w 



For radial motion v |l V</?jv the 3-acceleration is 



(| )FOT = _ (1 _3^ + 4^)V w 



2.3 PN equations of motion in GR 



In GR the equation of motion of a test particle is the geodesic equation [15] : 

^ = -r>V (20) 

where u 1 is 4-velocity of the test particle and rL is the Christoffel symbols. 

Post-Newtonian geodesic equation has been carefully studied in relativistic celestial 



mechanics and according to [21, 22] PN 3-acceleration of a test particle in the case of 



static spherically symmetric gravitational field is: 



:^)GR = -{1 + (1 + o)^ + (4 - 2a)^f - 3a(- • }V,:.v 



+ (4-2a)-(- • V<p N ) (21) 
c c 
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where v = dr/dt, tpjf = —GM/r, VtpN = GMr/r 3 . 

The very important quantity in equation (^l]) is the parameter of coordinate system 
a which depend on the choice of particular coordinate system and for example has the 
following values: a = 2 for Painleve coordinates, a = 1 for Schwarzschild coordinates, 
a = for harmonic or isotropic coordinates. Hence the trajectory of a test particle is a 
coordinate-dependent concept, i.e. it depends on the choice of a coordinate system due 
to parameter a. 

In the case of circular motion the 3- acceleration is: 

(^)gr = -{1 + (1 + a)^ + (4 - 2a)^-}V<p N (22) 
For radial motion the 3- acceleration is: 

~^)gr = -{1 - + (4 - 2a)— ^ 



(— )Gi? = -{l-3^ + (4-2a)^f}Vc^v (23) 



According to [21] in GR there are "coordinate-dependent (unmeasurable) quantities" 
and " measurable quantities" which may be directly obtainable from observation without 
involving any theoretical data (such as the laws of light propagation). For example the 
radius r of the circular orbit is a coordinate-dependent, unmeasurable quantity, while 
the masses M and m do not depend on coordinate system and hence are measurable 
quantities. 

The geodesic equations (^j] : 22, 23) are identical with equations (|l7| , p^JT9|) of FGT 
only for harmonic or isotropic coordinates where a = 0. Fock [23] also emphasized the 
special role of harmonic coordinates in GR. 



3 Trans lational motion of a rotating body in FGT 
3.1 PN equations of motion of a gyroscope in FGT 

Let us consider Poincare gravity force acting on a rotating body. ^From equation (p"2]) 
in the case of a gyroscope motion in static spherically symmetric gravitational field it 
follows the expression for the elementary gravity force dFp acting on each elementary 
mass dm of the gyroscope: 

dF P = -{(1 + ~ + 4^)V<^v - 3-(- • V<p N )}dm (24) 
2 c z c z c c 
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For a rotating solid body the total gravity force is the sum of elementary forces acting 
on elementary masses: 

F P = [ dF P (25) 



Taking into account that the velocity v of an element dm may be presented in the form 

v = V + [ur\ (26) 

where V is the translational velocity of the body, Q is the angular velocity, r is the radius 
vector of an element dm relative to its center of inertia, so that / rdm = 0. 
Inserting (|26"1) into (p4j) and (|25|) we get 



^_ M{(1 + f£ + ^ + f-g)V w 
V ,V -s 



_ 3 _(_ . v^) - — J [ur\([ur\ ■ V^ N )dm} (27) 

where M = J dm is the total mass of the body, I is its moment of inertia. 

Under the gravity force ( p7[ ) the rotating body will get the 3-acceleration according 
to general relation [15| : 



dv_jl-v*/<? {p _v { v p)) 



dt mo c c 
which may be written in the form: 

dV ^ V 2 ^ N 3/a; 2 , 

V ,V - , 3 



+4-(- • V VN ) + j^J [itf\<]pA ■ V<PN)dm (29) 



Equation of motion of a rotating body (|29| ) shows that the orbital velocity of the center 
of mass of the body will have additional perturbations due to rotation of the body. The 



last term in (29) depends on direction and value of the angular velocity u. It has an order 
of magnitude v 2 ot /c 2 and may be measured in laboratory experiments and astrophysical 
observations. 

3.2 Laboratory experiments with rotating bodies 



The most straightforward application of equation fl29| ) is to perform " Galileo-2000" ex- 
periment (which is an improved version of famous Stevinus-Grotius-Galileo experiment 



with free falling bodies in the Earth gravity field) just taking into account rotation of 
the bodies. 

Indeed let us consider three balls on the top of a tower (like the 110-m Drop Tower 
of the Bremen University) . The first ball is non-rotating and according to equation (p9|) 
its free fall acceleration is: 



m -(^)i = -(i-3^ + 4^)fy* (3U) 



Let the rotation axis of the second ball be parallel to the gravity force, i.e. u\\V<pn, 
hence its free fall acceleration is: 



3 R 2 uj 2 . 

92 = m + ^—) (3i) 



where it is taken into account that for homogeneous ball with radius R and mass M the 
moment of inertia is / = %MR 2 . 

5 

Let the rotation axis of the third ball be orthogonal to the gravity force, i.e. Co _L 
V(£>jv, hence its free fall acceleration is: 

^From equations (|30|),(|3l|),(^) it follows that considered three balls will reach the 
ground at different moments. True the difference is very small, for example if the radius 
of the ball is R = 10 cm and its angular velocity u = 10 3 rad/sec, then the expected 
difference in falling time from 110 m tower will be At « (l/2)(Ag/g)t ~ 2.5 x 10 -13 sec. 

For the NASA's gravity probe B (GPB) experiment |l8| the radius of the gyroscope 
is R ~ 2 cm and it rotates with about 2000 revolutions per second. The expected effect 
of the periodical perturbation of the acceleration of the gyroscope center of mass relative 
to non-rotating satellite center of mass is about Sg/g ~ 6 x 10~ 13 which is in principle 



measurable [17]. 



Another type of laboratory experiment for direct testing the velocity dependence of 



the Poincare gravity force (27) is to weigh a rotating bodies. If two bodies are at the 
balance and at a moment they start to rotate with different orientations of the rotation 
axes then the balance will be violated and hence measured by a scale. The expected 
difference in forces is again about uj 2 R 2 /c 2 . 
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In these laboratory experiments there are no problem with choosing a coordinate 
system at all. The height of a tower and the moments of the contact of the rotating 
bodies with the ground, and also readings of a balance scales are directly measurable 
quantities. Equation of motion ( |17| ) in field gravity theory gives the uniquely defined 
value for these laboratory experiments, while geodesic equation (^) includes arbitrary 
parameter a of a coordinate system. 

3.3 Periodical modulation of planet's orbits 

Post-Newtonian celestial mechanics in FGT does not depend on a choice of coordinate 
systems and is identical with PN approximation of GR only in harmonic or isotropic 
coordinates. This is why all classical post-Newtonian relativistic gravity effects are 



the same in FGT and GR. Note that according to [24] "astronomers have adopted the 
fairly standard convention of using "isotropic coordinates" rather than " Schwarzschild 
coordinates" when analyzing the solar system" (J24||,p.l097). 



According to [21] in GR the effects which depend on the choice of coordinate system 
(e.g. the parameter a) are unmeasurable. However there are such effects which are 
coordinate dependent in GR but are coordinate independent in FGT. An example of 
such effect is additional perturbation terms in equations of motion of rotating bodies 
which are determined by angular velocities of the bodies (see equation (|29|)). 

For a planet rotating with the angular velocity tJ in orbit with semi-malor axis a, and 
orbital angular momentum L the effect of rotation leads to periodical modulation of the 
orbital radius. The amplitude of the radius perturbation Aa may be roughly estimated 
as: 

Aa « a^- sin 2 9 (33) 
cr 

where R is the radius of planet, 8 is the angle between u) and L. 

According to equation ( |33"| ) the amplitude of the orbital modulation is about 6 cm for 
the Earth, 1 cm for Mars, 4.2 m for Jupiter and 340 m for Saturn. Future high accuracy 
radar observations of the solar system will test this tiny effects. 
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3.4 Rapidly rotating stars and pulsars in binary systems 

In a binary stellar system the rapidly rotating companion will produce the relativistic 
effect of the periodical modulation of the orbital motion. However for usual stars there 
are many other non-relativistic effects which usually hide small relativistic effects. 

The best candidates for testing relativistic gravity effects are pulsars in binary sys- 
tems. Timing observations of pulsars, i.e. the measurements of the time of arrival 
of pulsar signals at a radio telescope, are the high-precision experiments in astronomy. 
About 10% of the known pulsars are members of binary star systems, i.e. in orbit around 
a white dwarf, neutron star, or main-sequence star companion. Timing observations of 
these binary pulsars may be used for testing the periodical modulation of the orbit of 
the rotating bodies. 

In the case of equal masses of stars it is needed to generalize the above test-body 
approach to the case of two-body problem. However for rough estimation of the ex- 



pected effect it is enough to use equation (^). According to (33) the amplitude of the 
modulation of the orbit are determined by the spin angular velocity lu, the radius R of 
the pulsar (or more exactly, the distribution of the mass and velocity in the pulsar, i.e. 
its moment of inertia and spin angular momentum S), semi-major axis a and the 
angle p between pulsar's spin and direction of orbital angular momentum. 

Because the amplitude of the modulation effect is proportional to the semi-major 
axis, the largest effect may be observed for the special class of the long-orbital period 
binary pulsars |25|| . For known lj = 2tt/P p , a and R it is possible to determine the 
angle 9 P . If for a pulsar there is sufficiently accurate estimation of the angle 8 P , then it 
is possible to estimate the radius of the pulsar R. 

Let us make estimation for the expected effect in the case of two well-studied binary 
pulsars PSR B1913+16 and PSR B1259-63. The amplitude of the orbital modulation is 
given by (^3|) and the value of the expected variation in arrival times of pulses At = Aa/c 
may be written in the form: 



At « 1.75 • 10- 5 ( S ec)(f )(^) 2 ( T -^) 2 sin^ 9 P (34) 
Is F p 10 km 

where x = a ■ sin i/c is the projected semi-major axis, i is the inclination angle of the 
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orbit, P p is the period of the pulsar, R is the radius of the neutron star, and 6 P is the 
misalignment angle of the pulsar spin to the orbital momentum vector. 

For PSR B1913+16: period P p = 59 ms, x = 2.34 s, and 9 P « 22° @, hence the 
expected periodical variations in arrival time due to rotational effect is about 6 fxs. For 
PSR B1259-63 we have P p = 47.8 ms, x = 1296 s §|] and for = w ° the time variations 
will be about 0.7 ms. These numbers are close to the time residuals existing for these 
pulsars and it shows possibility to test the discussed effect. 

4 Conclusions 

Using Lagrangian formalism of relativistic field theory it is derived PN equations of 
motion of a rotating body in the framework of Feynman's field gravity theory. 

The trajectory of a test particle does not depend on a choice of a coordinate system 
and this is an essential property of field gravity theory which allow to make uniquely 
defined predictions for terrestrial laboratory experiments with rotating masses in Earth's 
gravity field. In these laboratory experiments there are no problem with choosing a 
coordinate system at all. The height of a tower and the moments of the contact of 
the rotating bodies with the ground, and also readings of a balance scales are directly 
measurable quantities. Equation of motion (|i~7| ) in field gravity theory gives the uniquely 
defined value for these laboratory experiments, while geodesic equation (]21j) includes 
arbitrary parameter a of a coordinate system. 

In the case of astronomical observations of distant objects it is important to consider 
also effects of emission, propagation and detection of photons, however it does not change 
the concept of coordinate independence of the trajectory of a test particle. 

Post-Newtonian periodical modulation of planetary motion in the solar system may 
be measured by future radar observations. Timing observations of pulsars in binary 
systems is the best test of this relativistic gravity effect. The expected arrival time 
variation due to pulsar rotation lies in interval from 1 /is up to 1 ms depending mainly 
on the value of the semi-major axis of pulsar orbit, misalingnment angle of pulsar spin, 
angular velocity of its rotation and the radius of a neutron star. In the case of binary 
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pulsars with known spin orientations this effect gives a possibility to measure radiuses 
of neutron stars. 
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